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Robert  Pinkel 


University  of  California,  San  Diego 
Marine  Physical  Laboratory  of  the 
Scripps  Institution  of  Oceanography 
San  Diego,  California  92132 


ABSTRACT 


Profiles  of  temperature  versus  depth  in  the  top  440  m of  the 
sea  were  taken  repeatedly  at  three  horizontal  locations  surrounding 
the  Research  Platform  FLIP.  The  time  fluctuation  of  the  temperature 

profiles  was  used  to  determine  isotherm  displacement  and  slope  varia- 
tion in  the  region  60-400  m. 

Measurements  were  taken  during  three  FLIP  operations.  Two 
were  off  the  California  coast,  in  November  1972  and  June  1973.  The 
third  was  north-northwest  of  Hawaii  in  November  1973. 

Horizontal  and  vertical  coherence  measurements  indicate  that 
the  internal  wavefield  can  be  divided  into  two  frequency  regions. 
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Above  2 eph  the  bandwidth  of  energetic  horizontal  and 
vertical  wavenumbers  is  comparatively  narrow.  The  isotherm  dis- 
placement spectrum  and  slope  spectrum  have  irregular  slopes.  A 
comparison  of  the  two  spectra  indicate  that  first  mode  is  strongly 
dominant. 

Below  2 eph,  the  bandwidth  of  energetic  wavenumbers  is 

broader,  corresponding  to  many  energetic  modes.  The  displacement 

_2 

spectrum  has  an  ft  frequency  dependence  down  to  about  two  cycles 
per  day.  The  slope  spectrum  has  an- approximate  spectral  form 

in  this  region.  This  indicates  that  the  bandwidth  of  energetic  hori- 
zontal wavenumbers  decreases  with  decreasing  frequency,  corresponding 
to  a bandwidth  of  vertical  wavenumber  approximately  constant  with 
frequency.  A vertical  spectral  analysis  of  the  data  substantiates 
this  general  picture,  while  revealing  the  presence  of  a low  frequency 
high  mode  region  of  surprising  spectral  intensity.  This  "shoulder" 
region  might  be  the  spectral  signature  of  horizontally  advecting 
features  in  the  temperature  field  which  are  not  necessarily  associat- 
ed with  vertical  density  perturbations.  However,  estimates  of 
horizontal  wavenumber  indicate  that  these  features,  as  well  as  those 
of  longer  wavelength,  tend  to  satisfy  the  linear  internal  wave 
dispersion  relation. 
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CHAPTER  I 
INTRODUCTION 

The  presence  of  internal  wave-like  motions  in  the  deep  sea 
has  long  been  known.  The  role  these  motions  play  in  the  energy  trans- 
fer and  mixing  of  the  ocean  is  still  not  clear. 

The  linear  mathematical  theory  of  internal  wave  propagation 
was  first  presented  by  Stokes  in  1847  and  has  been  thoroughly  explored 
since  then.  Extensive  theoretical  studies  of  the  interaction  of 
internal  waves  with  themselves  and  the  ocean  environment  have  been 
made.  Additional  understanding  awaits  meaningful  oceanic  measure- 
ments. 

Observational  efforts  have  yet  to  provide  a clear  picture  of 
the  internal  wave  field  in  the  open  sea.  The  wide  range  of  special 
and  temporal  scales  of  internal  waves  make  a single,  all  inclusive 
measurement  effort  impractical.  The  experiments  that  have  been  per- 
formed have  generally  been  of  two  classes.  Very  high  quality  single 
instrument  experiments  have  been  conducted,  either  by  towing  an  in- 
strument to  obtain  a spatial  data  series,  or  letting  it  drift,  to 
obtain  a time  series.  These  can  give  knowledge  of  the  amount  of 
motion  present  and  its  distribution  with  horizontal  or  temporal  scale, 
but  offer  only  indirect  indication  as  to  whether  internal  wave  dynam- 
ics truly  governs  the  motions.  A comparison  oc  the  temporal  and 
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three  dimensional  spatial  variability  of  the  notions  is  needed  to 
establish  the  relevance  of  the  linear  internal  wave  dispersion  rela- 
tion. Large  inul tidimensicnal  array  experiments  have  been  performed. 
Records  taken  by  different  instruments  in  these  arrays  have  generally 
shown  surprisingly  little  resemblance,  due  either  to  the  complexity 
of  the  wave  field  or  to  experimental  difficulties. 

The  measurements  to  be  discussed  are  from  a moderate  sized 
four  dimensional  (space-time)  array,  mounted  on  the  Research  Platform 
FLIP.  The  array  was  constructed  as  large  0s  practicable  while  still 
being  able  to  deploy  very  low  noise  sensors.  With  this  four  dimen- 
sional coverage  it  is  possible  to  get  a detailed  view  of  the  motions 
over  a small  area  of  the  ocean  around  FLIP,  with  consequent  poorer 
resolution  of  the  larger  scale  motions.  This  thesis  will  discuss  the 
results  of  the  first  series  of  operations  of  this  experiment.  The 
objectives  of  this  preliminary  study  are  to  categorize  the  distribu- 
tion of  energy  with  space  and  time  scales  and  to  ascertain  the  appli- 
cability of  linear  internal  wave  theory  to  the  motions  observed  in  the 
upper  ocean. 
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CHAPTER  II 

DESCRIPTION  OF  EXPERIMENT 


Experimental  Apparatus 

This  experiment  utilizes  temperature  as  a tracer  to  detect 
vertical  displacements  associated  with  internal  waves.  The  basic 
data  set  consists  of  profiles  of  temperature  versus  depth  in  the  upper 
440  m of  the  sea,  sampled  repeatedly  at  three  horizontal  locations 
surrounding  FLIP. 

The  primary  sensor  used  was  the  F-M  Thermometer  developed  by 
Mr.  Earl  Squier  (1967).  The  probe  consisted  of  three  thermistor 
beads  incorporated  into  the  feedback  loop  of  a phase  shift  oscillator 
circuit.  The  changing  resistance  of  the  beads  resulted  in  a varying 
output  frequency.  These  thermometer  oackages  were  characterized  by 
their  high  sensitivity  (~.003  °C  r.m.s.  typically)  and  minimal  low 
frequency  drift  ( ~ .05  °C  r.m.s.  over  several  weeks  of  operation). 

The  "rise  time"  of  the  thermometers  was  .25  sec. 

Each  of  the  3 profiling  units  consisted  of  several  components. 
An  F-M  thermometer  and  a vibratron  pressure  sensor  v/ere  mounted  on  a 
streamlined,  weighted,  five  feet  long  stainless  steel  probe.  The 
probe  was  lowered  into  the  sea  by  3/16"  double  armored  steel  cable. 

One  hundred  and  ninety  meters  up  the  cable,  a second  F-M  thermometer 
was  attached.  The  profiler  was  then  lowered  until  the  upper 
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thermometer  was  at  the  base  of  the  mixed  layer  (typically  65  m). 

From  this  position,  a remotely  controlled  five  horsepower  electric 
winch  cycled  the  thermometers  down  and  up  through  190  meters  at  a 
rate  of  4 m/sec.  This  yielded  a temperature  profile  of  the  top  380 
meters  of  the  thermocline  (Fig.  3).  Profiles  were  repeated  every  two 
minutes.  Measurements  were  taken  while  the  probes  descended.  Prelim- 
inary analysis  was  done  while  they  were  being  raised. 


To  produce  the  required  horizontal  separations,  it  was  necess- 
ary to  obtain  lightweight  booms  which  could  be  mounted  on  and 
deployed  from  FLIP.  Two  aluminum  girder  booms,  seventy-five  feet  long, 
were  constructed  and  installed  on  the  port  and  aft  faces  of  FLIP'S 
superstructure.  The  starboard  boom,  sixty  feet  in  length,  was  a con- 
verted aluminum  spinnaker  pole  from  the  schooner  GOODWILL.  When 
deployed,  a triangular  array  with  separations  38,  39  and  44  meters 
was  formed  (Figs.  1,  2).  The  size  of  the  array  was  limited  by  practi- 
cal considerations  associated  with  the  handling  of  the  booms  at  sea. 

The  azimuthal  positioning  of  the  array  was  affected  strongly 
by  the  wind,  which  tended  to  orient  FLIP  such  that  the  keel  was  20° 


to  the  left  of  the  wind  direction.  To  maintain  a constant  heading  for 
sufficient  time  to  gather  internal  wave  directional  information,  an 
automatic  orientation  control  system  was  employed.  This  consisted  of 
two  hydraulically  actuated  off-axis  propellers  which  were  mounted  on 
fLJP  s hull  and  controlled  by  the  ship's  gyro-compass.  Given  the 
consideiaole  viudage  of  the  array,  the  orientation  system  was  over- 
powered by  winds  over  ~10  knots.  However,  with  FLIP  in  the 
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Figure  1.  FLIP  internal  wave  array  configuration 


Pinkel 


equilibrium,  into-the-wind  position,  the  orientation  system  could 
reduce  the  azimuthal  fluctuations  to  approximately  + 8°.  Major  wind 
shifts  and  squalls  would  necessitate  the  repositioning  of  the  array. 

In  addition  to  the  basic  temperature  data,  supplemental  infor- 
mation was  recorded  from  a variety  of  sensors  which  were  borrowed 
from  others.  An  instrument  package,  incorporating  a sound  velocimeter, 
thermistor,  vibratron,  Doppler  flow  speed,  and  flow  direction  sensor 
was  loaned  by  Mr.  Earl  Squier  of  the  Marine  Physical  Laboratory  of 
Scripps  Institution  of  Oceanography  for  use  during  all  of  our  opera- 
tions. With  this  device,  necessary  vertical  profiles  of  horizontal 
currents  and  density  could  be  obtained.  Wind  speed,  direction  and 
atmospheric  pressure  were  recorded  during  portions  of  June  1973  and 
November  1973  FLIP  cruises  using  instruments  of  Dr.  Carl  Friehe 
(Department  of  AMES,  University  of  California,  San  Diego).  Surface 
wave  measurements  were  conducted  during  the  November  1972  and  June 
1973  trips  with  sensors  supplied  by  Dr.  Hal  Guthart  and  Mr.  Bob  Martin 
of  Stanford  Research  Institute. 

Data  handling  functions  were  performed  at  sea  by  a Hewlett 
Packard  2116B  computer  system  developed  by  Dr.  Bruce  Williams,  Marine 
Physical  Laboratory,  Scripps  Institution  of  Oceanography.  Fifteen 
F-M  and  five  analog  data  channels  were  scanned  while  the  probes  fell. 
Incoming  information  was  stored  in  memory  and  written  on  nine  track 
magnetic  tape  immediately  after  the  probes  completed  the i r descent. 

All  at-sea  analysis  and  data  display  tasks  were  performed  v/hilc  the 
probes  were  rising,  between  drops.  During  tin's  period  the  computer 
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performed  a series  of  interpolations,  picking  out  the  set  of  depths 
at  which  a pre-selected  set  of  temperatures  had  occurred  in  the  pro- 
file of  each  of  the  thermometers.  When  temperature  inversions  were 
present,  the  shallowest  depth  at  which  the  given  temperature  was 
encountered  was  arbitrarily  chosen.  These  isotherm  depths  were  then 
stored  in  memory  and  later  written  on  tape  along  with  the  profile  data 
from  the  next  drop.  Selected  isotherm  depths  from  each  probe  were  also 
stored  on  the  computer^  disk  for  display  and  analysis  purposes  at  sea. 
The  disk  held  approximately  24  hours  isotherm  data.  The  isotherm 
interpolation  required  about  ten  seconds  time  during  the  probes' 
ascent;  the  roughly  forty  seconds  time  remaining  was  used  for  display 
of  the  isotherm  information  or  raw  data  on  the  computer's  cathode  ray 
storage  scope. 

Data  Collection 

Measurements  were  made  during  three  FLIP  cruises.  The  first 
two  were  conducted  off  the  California  coast  during  November  1972  and 
June  1973.  The  starting  location  for  both  was  31°N,  121°W,  a region 
of  four  km  deep  water,  with  flat  bottom  and  scattered  seamounts.  This 
site  is  120  km  offshore  (southwest)  of  the  Patton  Escarpment,  a three 
km  step-like  break  in  the  continental  shelf.  In  this  region,  the 
southward  flowing  water  of  the. Cal ifornia  Current  intermix  with  the 
northward,  more  intermittent  flow  of  the  Davidson  Current. 

The  Novemb'  1972  operation  lasted  three  weeks.  It  was  in- 
terrupted briefly  when  (LIP  drifted  eastward  toward  the  continental 
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shelf  and  had  to  he  re-towed  into  deep  water.  On  the  June  1973  cruise 
FLIP  was  on  station  four  weeks.  The  weather  was  quite  varied  on  both 
cruises.  Maximum  winds  were  typically  25-35  knots.  Periods  of  calm 
were  frequent.  Pronounced  weather  changes  occurred  on  a four  day  - 
one  week  time  scale. 

It  was  not  known  how  representative  these  California  measure- 
ments were  of  the  wave  field  in  the  central  areas  of  the  deep  sea. 

The  proximity  of  the  continental  shelf  (5-15  internal  tidal  wave- 
lengths distant)  was  one  cause  of  concern.  The  shelf  has  been  ident- 
ified as  a possible  source  (Kattray,  et  al . , 1969;  Wunsch  and  Hendry, 
1972)  and  sink  (Wunsch,  1968)  of  internal  wave  energy.  The  complex 
horizontal  current  structure  of  the  region  could  also  be  expected  to 
influence  internal  wave  propagation.  As  a result,  it  was  decided  to 
make  a short  cruise  north/northwest  of  Hawaii  during  early  November 
1973.  Eleven  days  of  data  were  taken  starting  at  25°N,  159°W  in  four 
km  deep  water.  The  sea  floor  was  smooth  with  no  seamounts  in 
the  immediate  area.  Conditions  were  quite  moderate  for  the  first 
week,  and  rough  with  25  knot  winds  for  the  remaining  four  days.  In 
all,  approximately  fifty  days  of  usable  data  were  collected.  During 
that  time  at  lea.<t  one,  and  typically  all  three,  profilers  were  in  op- 
eration. Elec  rical  problems  witli  underwater  connectors  were  the  most 
common  difficulty  encountered.  Severe  squalls  occasionally  necessita- 
ted shutdown  of  the  aft  profiler.  During  those  high  winds  El  IP’s 
drift  would  become  sufficient  to  cause  the  instruments  to  be  swept 
into  the  hull.  Squalls  would  also  overpower  El  IP's  orientation 
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system,  azimutnally  spinning  the  ship  and  interrupting  time  series 
measurement  of  the  directional  properties  of  the  wave  field.  Power 
generation  difficulties  aboard  FLIP  and  at  sea  refueling  were  rare 
causes  of  complete  cessation  of  the  experiment. 

Profiles  of  the  Vaisala  Frequency 

Representative  profiles  of  the  Vaisala  frequency 

N(z)  - (2-  (|f-  - £2.))1/2 

p az  Cz 


for  the  cruises  are  presented  in  Figure  4.  Typically,  the  November 
1973  Hawaii  cruise  had  the  greatest  Vaisala  frequency;  the  June  1073 
California  cruise  the  smallest.  These  profiles  are  composites, 
resulting  from  many  sound  velocity-temperature-depth  lowerings  from 
FLIP.  An  STD  system  was  also  used  in  the  June  1973  cruise.  These 
profiles  will  be  utilised  to  tune  the  analysis  procedures  to  be 
described. 

It  was  found  that  the  Vaisala  profile  in  the  upper  thermocline 
changes  significantly  on  a time  scale  comparable  to  that  of  low  fre- 
quency internal  waves.  An  extreme  example  of  the  degree  of  modifica- 
tion in  the  upper  thermocline  which  can  occur  is  given  in  Figure  5. 
Thirteen  hours  of  temperature  profile  data,  taken  on  November  16,  1972 
are  presented.  Successive  profiles  are  separated  by  four  minutes. 

The  severe  modification  in  the  temperature  structure  above  120  meters 
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over  the  thriteen  hour  series  is  striking.  The  many  short  lived 

inversions  visible  suggest  that  vertical  mixing  processes  might  be 
active  in  this  region. 

The  migrating  zones  of  locally  high  temperature  gradient  also 
stand  out.  It  is  doubtful  that  these  featu-es  are  a result  of  small 
vertical  scale  internal  waves,  straining  an  otherwise  smooth  tempera- 
ture field.  The  profile  is  frequently  completely  deformed  into  a 
dynamically  stable  temperature  inversion,  which  cannot  possibly  result 
from  internal  wave  vertical  displacements. 

Mixing  and  advective  effects  are  clearly  modifying  these 
profiles.  It  must  be  assumed  that  their  influence  on  the  density 
field  will  be  Significant  also.  Small  scale  details  in  the  density 
profile  presumably  evolve  more  rapidly  than  the  large  scale  features. 
As  a consequence,  a vertical  length  scale  .an  be  defined  under  which 
the  dansity  profile  significantly  changes  over  the  time  period  necess- 
ary to  obtain  statistically  stable  measurements  of  the  wave  field. 

The  analysis  method  to  be  utilized  will  fail  to  characterize  motions 
of  vertical  scale  less  than  this  satisfactorily.  It  will  be  seen 
that  noise  in  the  analysis  resulting  *rom  the  time  variation  of  the 
Vaisala  profile  and  experimental  error  in  its  estimation  will  signifi- 
cantly contaminate  the  measurements  of  motions  of  less  than  40  meter 
vertical  wavelength. 

Poscri ptiun  of  the  Data 

A thirteen  hour  segment  of  isotherm  data  from  the  starboard 
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profiler  is  presented  in  Figure  6.  This  sample  is  not  representative 
of  the  bulk  of  the  observations.  However,  it  contains  examples  of 
most  of  the  features  found  in  the  data.  The  isotherms  followed  by 
each  thermometer  are  chosen  at  equal  increments  of  temperature  (~.  1°C 
for  the  upper  thermometer  and  ~.05CC  for  the  lower).  Close  spacing  of 
the  isotherms  indicates  a high  vertical  temperature  gradient.  Gaps 
represent  isothermal  regions  or  temperature  inversions.  The  large 
gap  from  220  to  280  meters  lies  between  the  lowest  isotherm  followed 
by  the  upper  thermometer  and  the  highest  isotherm  followed  by  the 
lower. 

The  generally  smooth  nature  of  the  time  oscillations  of  the 
isotherms  can  be  clearly  seen.  This  implies  a relative  absence  of 
energy  in  the  motions  at  frequencies  high  compared  to  the  Nyquist 
frequency  of  15  cph.  Exceptions  are  apparent  where  isotherms  jump 
frcm  top  to  bottom  of  nearly  isothermal  zones  and  visa  versa.  These 
jumps  are  due  to  slight  temperature  changes  within  the  zones  (perhaps 
associated  with  horizontal  advection).  They  clearly  do  not  represent 
vertical  motions  of  the  water  and  can  be  easily  corrected,  using 
vertical  velocity  information  from  the  isotherms  above  and  below  the 
jump. 

Considering  now  only  the  smoothly  oscillating  component  of 
the  data,  two  features  are  quite  apparent.  The  higher  frequency 
motions  appear  to  be  quite  persistent  with  depth.  Local  crests  present 
at  eighty  meters  can  frequently  be  seen,  with  the  same  phase,  one 
hundred  fifty  motors  below.  In  contrast,  the  lower  frequency 
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oscillations  tend  to  exhibit  greater  variation  with  depth.  Only  at 
frequencies  less  than  3 cph  are  large  phase  changes  in  the  vertical 
noticeable.  In  many  records,  the  M2  internal  tide  can  be  seen  to 
undergo  a vertical  phase  reversal  in  less  than  100  m.  It  is  clear 
from  this  record  that  vertical  wavelengths  long  compared  to  several 
hundred  meters  dominate  the  wave  field  at  frequencies  greater  than 
3 cph.  At  lower  frequencies,  shorter  vertical  wavelengths  are  associ- 
ated with  the  observed  vertical  variations. 

Fundamental  Experimental  Uncertainties 

The  isotherm  fluctuations  will  be  used  to  infer  changes  in 
the  r.m.s.  potential  energy 

Ep  = 1/2  p nhz)  <„2> 

3 

p = deniity  gm/cm 
N(z)  = Vaisala  frequency 

n = displacement  of  an  isopycnal 
surface  above  its  mean  position 

and  the  r.m.s.  vertical  kinetic  energy 
Ek  = V2  p<w2> 

w = vertical  velocity 

in  the  upper  thormocl ino.  Since  temperature  rather  than  density  is 
measured,  error  is  introduced  into  the  potential  energy  determination 
whenever  temperature  becomes  a poor  tracer  of  deniity.  Mixing, 
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diffusion,  and  lateral  advection  are  the  principal  sources  of  this 
error.  Uncertainty  in  the  vertical  kinetic  energy  estimation  is 
introduced  by  the  fictitious  vertical  velocity 

e(z»t)  = - uTzTtl  • yj/ll 

n 9Z 

which  results  from  the  horizontal  advection  of  the  water  through  the 
array.  This  effect  would  be  present  even  if  density  and  not  tempera- 
ture were  measured.  It  is  traditionally  discussed  in  terms  of  the 
Doppler  shifting  of  measured  frequencies. 

Neither  error  can  be  a priori  distinguished  from  vertical  in- 
ternal wave  velocities  using  temperature  data  alone.  Yet  the  role 
they  play  is  not  mysterious.  Their  effect  on  the  interpretation  of 
the  measurements  is  discussed  below. 

The  potential  energy  error  results  from  temperature  being  a 
poor  tracer  of  vertical  motions.  Figure  6 reveals  several  examples 
suggestive  of  mixing  or  horizontal  intrusive  events  causing  vertical 
fluctuation  of  the  isotherms.  Consider  the  motions  occurring  at 
depth  325  nr,  time  8-12  hours.  Isotherm  after  isotherm  is  seen  to 
migrate  from  top  to  bottom  of  a region  of  low  temperature  gradient, 
possibly  associated  with  the  advance  of  a warm  lens  of  water  into  the 
FLIP  sampling  area.  Note  that  all  motions  of  this  sort  in  Figure  6 
are  typically  confined  to  an  rsolated  10-25  m depth  zone.  The  event 
dominates  the  signal  at  one  depth  but  has  little  apparent  effect  on 
the  isotherm  motions  30  or  40  meters  away.  Thus,  given  the  350-400 
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meter  vertical  picture  from  the  array,  it  is  seldom  difficult  to  sep- 
arate out  the  longer  vertical  scale  wave  motions  from  the  local  intru- 
sive mixing  effects.  However,  there  is  no  way  to  separate  the 
shorter  (~20  m)  vertical  scale  wave  motions  from  the  more  gradual, 
less  "obvious"  intrusions  without  direct  density  measurements.  As  a 
result,  a detailed  study  of  the  short  vertical  scale  fluctuations  will 
not  be  attempted. 

Error  in  the  vertical  kinetic  energy  estimates  is  associated 
with  the  Doppler  shifting  of  measured  frequencies.  If  u is  the 
true,  intrinsic  frequency  of  a wave,  and  q the  observed  Doppler 
shifted  frequency,  the  relative  frequency  error  610  is  given  by 


£ _ Q— U1  _ 

ui 


u(z) 


where  Cp  is  the  horizontal  phase  velocity  of  the  wave.  As  internal 
wave  horizontal  phase  velocities  are  not  large  compared  to  typical 
relative  advective  velocities  (u(z)«10  cm/sec)  significant  Doppler 
shifting  might  be  expected. 

In  principle,  with  the  four  dimensional  data  set  provided  by 
this  experiment  and  the  assumption  of  an  internal  wave  dispersion 
relationship,  the  intrinsic  frequency  could  be  obtained  from  measure- 
ments of  the  horizontal  and  vertical  spatial  frequencies.  Success  in 
this  approach  is  curtailed  both  by  the  limited  resolution  of  the  array 
and  the  inapplicability  of  the  theoretical  linear  dispersion  relation, 
as  a result,  primarily  of  the  pronounced  current  shear  in  the  upper 
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thermocl ine.  The  data  will  thus  be  analyzed  and  displayed  in  terms 
of  the  encounter  frequency  n.  Measurements  taken  during  periods  of 
excessively  high  relative  FLIP  drift  will  be  discarded.  Where  signi- 
ficant wave  energy  is  found  associated  with  slower  phase  velocities 
it  must  be  horn  in  mind  that  it  has  possibly  been  severely  Doppler 
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CHAPTER  III 


THE  STATISTICAL  DESCRIPTION  OF  THE  MOTION 


The  isotherm  motions  measured  can  be  mathematically  modeled 


by 


2tt 


I (x ,z, t ) = 


dA(< , 0 ,u))W(k  ,oj , z)e^K  x 


K = 0 0“O  U)=0 


where  W is  a solution  to  the  linear,  shear  free  internal  wave 
equation 


+ / —■ W = 0 

az2  (a,  -r) 


W = 0 at  z = 0 


and  is  normalized  such  that 


_L 

2D 


p(z)  N2(z)  W2(K,u),z)dz  = 1 


This  representation  is  similar  to  the  equivalent  continuum 
approach  of  Garrett  and  flunk  (1972)  and  is  discussed  in  greater  detail 
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in  Appendix  I. 

The  r.m.s.  potential  energy  is  given  by 


Ep  = 1/2  „(z)N2(z)  <n2> 


2lT 

K=0  0 = 0 OJ  = 0 


w £.71  go 

= 1/2  p(z)N2(z)  j j j <dAdA*(K,0,w)>  W2(k,u,z) 


2tt 


oo  C7T  o° 


= 1/2  p (z)N2(z)  f ' ( dE(<,o,u)  W2(k,co,z) 

J 


K=0  0 = 0 Ci}-0 


where 

dE(<  ,o  ,w)  = <dAdA*(K  ,0  ,u)  )> 

is  the  potential  energy  spectrum.  The  objective  of  this  section  is 
the  estimation  of  the  spectrum. 

The  spectrum  can  be  related  to  the  data  set  through  the 
covariance  function 

(^(x-j-x^.Zi  = -.ri(xi  ,z-|  ,t-|  ^(x^.z^.t^ 

or  more  convert  l iently,  its  time  Fourier  transform,  the  cross  spectrum 
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CR(x^-X2>z^  - F(C(x^ -X2 ,z-|  ,Z2,At-|2)  ;At-| 2) 

<F (n (xi  >zi  >t] ) Jt'j ) F*(n  (X2,Z2  *’ ^2^ 


where  F(a;y)  represents  the  Fourier  transform  of  a with  respect 
to  y.  Ideally,  the  spectrum  can  then  be  obtained  through  Fourier 
transform  of  the  horizontal  dependence  of  the  cross  spectrum 

A A 

W(k,w,z-|)  W(k,w,Z2)  dE(<,co,o)  = F(CR(ax,z-|  ,Z2»uj) ; >0  . 


Alternatively,  a wavenumber  magnitude-frequency  spectrum  results  from 
vertical  integral  transformation  of  the  cross  spectrum,  utilizing  the 
internal  wave  vertical  wavefunctions  introduced  above. 


A 

dE(tc  ,00 ) 


2n  -H 

dE(ic,e,w)  = j dz-| 
0 0 


/N^(z1)-w 


CO 


2 


-H  w2/  v 2 
N (z2)-co 

W(k,u,Zi  ) dZ2 ( 1 2'  1 2 ^ W(k,co,Z2)  CR(o,z-|,Z2,co)  . 

co  -f 

0 


* The  equivalence  of  these  two  approaches  depends  on  the  assump- 
tion that  linear,  snear  free  internal  wave  dynamics  well  relates  the 
horizontal,  vertical  and  temporal  scales  of  the  observed  motions. 

Given  the  typical  magnitude  of  horizontal  shears  in  the  thermocline, 
the  validity  of  this  assumption  is  doubtful.  An  effort  will  be  made 


25 


--->•■  — ■ . - • ~ ~ .... ■■■■  - A, 


Pinkel 


to  compare  these  approaches  at  the  end  of  this  chapter. 

This  experiment  does  not  provide  sufficient  information  to 
perform  either  of  the  above  integral  transform  direct  spectral  esti- 
mation methods.  Nevertheless,  a sample  cross  spectrum 
CR(x.j -Xj ,z.j  ,Zj  ,n ) can  be  evaluated  from  the  data  obtained,  and  less 
direct  methods  used  to  infer  the  form  of  dE(K,o,w).  The  various 
spectral  estimators  to  be  used  in  this  study  will  all  be  related  to 
CR(x\ ,zi  ,z  . ,n). 

To  calculate  the  sample  cross  spectral  estimate,  the  measured 
region  of  the  water  column  was  divided  into  forty  fixed  depth  inter- 
vals. For  every  drop,  the  mean  internal  wave  vertical  displacement 
was  calculated  in  each  of  the  depth  intervals  at  each  of  the  three 
piofileis.  Forty  x 3 time  series  of  mean  isotherm  displacement  were 
thus  produced.  These  fixed  depth  series  were  divided  into  17.07  hour 
(512  drop)  records,  first  differenced,  multiplied  by  a triangular  data 
window,  and  fast  Fourier  transformed.  The  resulting  spectral  window, 
sine  ( 2a f / , had  elemental  bandwidth  2/17  cph.  Twenty-three  overlapping 
17  hour  series  were  analyzed  from  the  November  1972  cruise,  thirty  for 
the  June  1973  cruise  and  seventeen  from  the  November  1973  operation. 

In  the  sections  that  follow,  the  sample  spectral  estimators 
will  be  presented  and  their  relation  to  the  spectral  density  function 
dE(K,o,ui)  discussed.  First  the  frequency  dependence  will  be  examined, 
then  the  frequency  and  depth  dependence;  then  frequency  and  horizontal 
dependence.  Finally,  the  four  dimensional  behavior  of  the  spectrum 
will  be  investigated. 
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T_he  Vorti_caj_  Pi  sjl  acement  Frequency  Fppr t rum 

The  vertical  displacement  frequency  spectrum 


= CR(o,z,z,a) 


m2/cph 


is  related  to  the  spectral  density  function  by 

» 2tt 

^2,“^  [ I dE(K,Q>aj)  W2(K,^,z). 

K=0  0=0 


For  a given  frequency  and  depth,  the  value  of  the  displacement 
spectrum  is  given  by  the  ~ - integrated  sum  of  all  energy  at  that 
frequency,  modulated  by  the  square  of  the  vertical  wavefunctions 

W (r:*Ji,z)  evaluaLed  at  the  given  depth.  For  frequencies  far  above 
the  inertial  frequency,  S(z,fl)  is  directly  proportional  to  the 
potential  energy  spectrum  (energy/vol/cph)  with  proportionality 
constant  1/2  p(z)N2(z).  The  spectra  to  be  presented  below  are  thus  an 

indicator  o.  the  distribution  of  internal  wave  potential  energy  with 
frequency. 

A displacement  spectrum  of  the  20°  isotherm  taken  from  10  days 
of  the  Hawaii  data,  is  presented  in  Figure  7.  Below  2 cph  the  spec- 
trum can  be  approximately  characterized  by  an  jf2  slope.  Between 
2 q>h  a'Kl  thc  1ocal  '^isala  frequency  there  is  frequently  some  irr-  gu- 
1 n r i ty  in  the  slope;  in  this  case  a peak  is  present.  At  the  Vaisala 
frequency,  the  spectrum  cuts  off  sharply,  falling  15  to  25  dB  in 


Frequency,  c.p.h. 


Figure  7.  Hawaii  isotherm  displacement  spectrum,  40  degrees  of  freedom 
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energy  down  to  a spectral  region  dominated  by  thermometer  noise.  The 
slope  of  the  spectral  cutoff  is  typically  rf10  (however  the  spectral 
processing  techniques  used  were  not  intended  to  produce  an  optimal 
spectral  estimate  in  the  region  of  the  cutoff.  The  sf10  slope 
must  be  considered  an  imprecise  estimate).  The  high  frequency  noise 
level,  if  assumed  to  be  white  noise,  corresponds  to  an  r.m.s.  uncer- 
tainty in  the  depth  of  an  isotherm  of~30  cm.  Lower  noise  measure- 
ments taken  in  June  off  California  indicate  a r.m.s.  precision  of 
~15  cm,  which  approaches  the  point  where  the  concept  of  a smoothly 
varying  isothermal  surface  is  no  longer  viable  in  the  upper  ocean. 

The  absence  of  a pronounced  spectral  cutoff  in  many  previous 

internal  wave  experiments  utilizing  fixed  depth  instruments  has  been 

attributed  by  Phi  1 1 i ps  (1971)  and  Garrett  and  Munk  (1971)  to  spectral 

contamination  effects  associated  with  nonlinearities  in  the  vertical 

gradients  of  the  properties  measured.  Assuming  a step-like  vertical 
-2 

gradient,  an  sj  contamination  spectral  form  was  indicated  by  their 
studies.  The  level  of  contamination  could  be  sufficient  to  dominate 
the  signal  in  the  internal  wave  band,  in  addition  to  obscuring  the 
cutoff 

To  verify  their  suggestion,  a sub-experiment  was  conducted 
during  a July  1971  cruise  off  the  California  ccast.  A temperature 
sensor  was  lowered  to  a fixed  depth  (80  m)  while  a profiler  was  in 
operation  10  meters  away.  Eight  hours  of  data  were  taken  to  provide 
a 16  degree  of  freedom  look  at  the  high  frequency  end  of  tiie  spectrum. 
The  displacement  spectrum  of  the  isotherm,  figure  8,  shows  a 
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H.  California  isotherm  displacement  spectrum.  July  1971, 
record,  ~45  cph  in' yquist  frequency,  16df. 
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pronounced  peak  prior  to  a severe  ~40  dB  Vaisala  cutoff.  The  spectrum 
of  temperature  fluctuations.  Figure  9,  displays  a noticeable  ~10  dB 
Vaisala  cutoff,  followed,  at  higher  frequency,  by  the  predicted  n 
contamination  spectrum.  The  "fine  structure  contamination"  did  not 

% 

dominate  the  signal  in  the  internal  wave  frequency  band  during  this 
observation.  Recently  Briscoe  (personal  communication)  has  obtained 
fixed  depth  temperature  spectra  from  a TRI -mooring  in  the  Atlantic 
Ocean  which  also  appear  to  be  relatively  uncontaminated  by  nonlineari- 
ties in  the  temperature  gradient  and  are  reminiscent  of  the  isotherm 
displacement  spectrum  presented  above. 

The  spectral  peak  prior  to  the  Vaisala  cutoff  has  been  a 
source  of  great  puzzlement.  An  explanation  has  been  suggested  by 
Desaobies  (pers.  comm.)  and  Garrett  and  flunk  (pers.  comm.)  which 
can  be  loosely  paraphrased  as  follows:  Assuming  a Vaisala  profile 

monotonically  decreasing  with  depth,  for  a given  frequency,  u,  one 
can  define  a cutoff  depth,  D(w)  such  that  N(D(00) ) = u.  Below  this 
depth,  wave  functions  W(K,u,z)  decay  exponentially  for  all  *. 
Slightly  shallower  than  this  "turning  depth",  depending  on  K,  all  wave 
functions  will  have  a relative  maxima,  exhibiting  a behavior  similar  to 
the  Airy  function.  At  depths  further  into  the  wave  guide,  the  phases 
of  the  vertical  functions  will  tend  to  de-correlate,  due  to  the  dif- 
fering vertical  "wavelengths"  associated  with  differing  Recalling 
the  del  ini  Lion  of  the  displacement  spectrum 
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CD  2 71 

S(z.n)  = J | <dAdA*U,e,u)>  W2(K,n,z) 

o o 

< e 

we  see  that  the  < integration  will  tend  to  weight  some  regions  of 
wavenumber-frequency  space  more  than  others,  depending  on  whether  the 
corresponding  W (k.sz.z)  is  near  a node  or  a maximum.  However,  at 
depths  just  above  D(w),  all  wave  functions  should  be  near  a relative 
maximum.  The  resulting  constructive  interference  would  produce  a peak 
in  the  displacement  spectrum.  Correspondingly,  at  a given  depth,  the 
spectral  peak  should  occur  at  frequencies  slightly  below  the  Vaisala 
cutoff. 

This  simplified  discussion  rests  heavily  on  the  assumption 
that  <dAdA*(<,o,n)>  varies  smoothly  with  u and  *.  Otherwise,  it 
would  be  impossible  to  distinguish  between  the  effects  of  the  above 
mentioned  interference  phenomena  and  the  n,*  variation  of 
<dAdA*(K,0,n)>  itself. 

Figure  10  displays  the  same  displacement  spectrum  as  Figure  7, 
plotted  together  with  the  spectrum  of  the  15°  isotherm  (mean  depth 
380  m).  As  the  data  were  taken  at  the  same  time,  spectral  differences 
presumably  result  from  the  depth  variability  of  the  square  of  the 
vertical  wavefunction,  W (,;,q,z).  Several  differences  are  clearly 
apparent  . The  spectrum  of  the  colder  isotherm  tends  to  have  greater 
variance  than  tluil  ot  the  warmer  at  frequencies  below  the  Vaisala  cut- 
off. The  cu tot f occurs  at  a lower  frequency,  indicative  of  the  de- 
crease in  the  Vaisala  frequency  with  depth.  Also,  the  pronounced 
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spectral  peak  observed  at  90  m is  not  present  at  380  m. 

To  extend  this  comparison,  six  Hawaii  displacement  spectra 
(72  degrees  of  freedom)  are  plotted  in  Figure  11.  The  uppermost 
spectrum  is  plotted  with  the  correct  scale,  the  others  are  succes- 
sively shifted  by  10  dB  to  prevent  the  plots  from  overlapping.  The 
depth  variations  in  spectral  slope,  cutoff  slope  and  high  frequency 
peak  are  here  clearly  apparent.  Presented  in  the  same  format  are  dis- 
placement spectra  from  the  November  1972  California  (Fig.  12)  and  June 
1973  California  (Fig.  13)  cruises.  The  November  1972  spectra,  estima- 
ted at  64  degrees  of  freedom  are  the  least  energetic,  and  tend  to  have 
the  weakest  cutoffs.  The  pre-cutoff  multiple  spectral  peaks  which 
appear  to  be  statistically  significant  in  the  upper  spectra  "fade 
away"  by  the  lower  spectrum. 

The  more  energetic  June  California  data  (48  df)  show  no 
clearly  defined  pre-cutoff  peaks;  rather  a broad  rise  is  present 
(upper  spectra)  or  nothing  at  all  (lowermost  spectrum).  The  cutoff 
slopes  are  quite  pronounced,  and  the  high  frequency  noise  levels  are 
the  lowest  of  any  of  the  three  curises. 

Figures  14  and  15  are  overlays  of  the  preceeding  plots, 
which  permit  direct  comparison  of  the  late  fall/late  spring  California 
spectra  and  late  fall  California  and  Hawaii  spectra.  Note  that  in  the 
region  of  the  water  column  measured,  the  November  Hawaii  wave  field 
has  approximately  1.5  times  the  potential  energy  of  the  November 
California  wave  field. 


/cph 


Figure  12.  California  November  1972:  Isotherm  displacement  spectra  64  df. 
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Cycles/Hour 


Htfure  15.  California  - Hawaii  spectral  comparison. 
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The  Vertical  Velocity  Spectrum 

The  spectrum  of  isotherm  vertical  velocity  is  related  to  the 

2 

isotherm  displacement  spectrum  by  the  factor  ft 
V(z,ft)  = <|“  (x,z,ft)|2>  = ft2  S(z,ft)  . 

Spectra  of  vertical  velocity  for  the  three  cruises  are  pre- 
sented in  Figures  16-18.  Table  I lists  the  associated  variances.  The 

pre-cutoff  spectral  rise  in  the  displacement  spectrum  is  transformed, 

2 

by  the  ft  weighting,  into  a distinct  peak  in  the  vertical  velocity 
spectrum.  Frequently  half  of  the  total  vertical  velocity  variance 
measured  is  located  in  this  2-3  cph  wide  pre-cutoff  peak. 

At  greater  depths  this  velocity  spectral  peak  might  well  be 
more  localized  in  frequency.  It  will  be  suggested  that  more  inodes  are 
present  at  lower  frequencies.  Given  the  general  reduction  in  the 
Vaisala  frequency  with  depth,  more  modes  will  be  present  near  the 
cutoff  frequency  to  sharpen  the  spectral  peak.  Deeper  measurements  by 
Voorhis (1968) , Cairns  (in  press)  and  Briscoe  (pers.  comm.)  appear  to 
substantiate  this  assertion. 

Given  the  large  number  of  observations,  it  seems  reasonable  to 
conclude  that  internal  wave  induced  Vaisala  "ringing"  of  the  deep  sea 
is  a common  phenomena,  which  is  associated  with  up  to  half  the  verti- 
cal kinetic  energy  present  in  the  sea. 


(cm/sec  )2/cph 
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TABLE  I 

VERTICAL  VELOCITY  VARIANCE 


November 

■ 1972 

June 

1973 

November 

■ 1973 

Depth 

Variance 

Depth 

Variance 

Depth 

Variance 

m 

2 2 
cm  /sec 

m 

2 2 
cm  /sec 

m 

2 2 
cm  /sec 

88 

.2924 

104 

.2713 

94 

.7860 

118 

.3748 

140 

.3460 

134 

.8211 

158 

.3920 

182 

.3901 

182 

.6984 

255 

.2847 

219 

.3687 

276 

.5383 

298 

.3018 

329 

.3657 

332 

. 5301 

239 

.2663 

394 

.3498 

390 

.5808 
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The  Vertical  Cross  Spectrum 

The  depth  variability  of  the  wave  field  can  be  further  studied 
through  estimates  of  the  vertical  cross  spectrum 


V C(z1 , z2,n)  = CR(0,z-j  ,z2,n)  . 


Cross  spectra  are,  in  general,  complex  functions  of  frequency 
and  spatial  separation.  Considering  sinusoidal  waveforms,  Munk  et  al . 
(1964)  has  related  the  real  part  of. the  cross  spectrum  (co-spectrum) 
to  the  sum  of  the  energy  propagating  from  measurements  points  x^  to 
x2  and  x^  to  x-j . The  imaginary  part  (quadrature  spectrum)  was  re- 
lated to  the  difference  in  energy  propagation.  It  can  be  verified  that 
the  same  result  applies  for  internal  wave  vertical  functions.  It 
follows  that  our  wave  field  model,  which  assumes  no  net  vertical 
energy  flux,  requires  that  the  vertical  cross  spectrum  will  be  real 


VC(zi  ,z^,n)  ~ <r|(*o’^l 


n) n* ( ^ ) > 


2,i 

’ 

4 

0 

0 


:dAdA*(K.e,n)>  W(K,n,z1  )W(K ,r>,z2)  . 


Ttie  ratio  of  quadrature  spectrum  estimate  to  the  co-spectrum  estimate 
will  thus  serve  to  quantify  the  extent  to  which  the  assumption  of 
standing  waves  in  the  vertical  is  appropriate. 

To  increase  the  statistical  stability  of  the  vertical  cross 
spectral  estimate,  a depth  lag  averaging  technique  was  utilized.  In 
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a given  frequency  band,  averages  were  formed  of  cross  spectral  estima- 
tors of  like  value  of  the  W.K.B.  stretched  separation  depth 


producing  a final  estimate  which  is  a function  not  of  z^  and  z ^ 
but  only  the  stretched  separation  distance  AZ.  Implicit  in  this 
procedure  is  the  assumption  that,  aside  from  the  scale  distortion 
caused  by  the  variation  of  N in  the  vertical,  the  wave  field  is 
otherwise  statistically  stationary  with  depth. 

Contour  plots  of  the  vertical  co-spectrum  and  quadrature 
spectrum  of  the  Hawaii  data  are  presented  in  Figure  19.  Both  spectra 
are  normalized  by  Re(VC(o,n)).  The  quadrature  spectrum  varies  irreg- 
ularly, rarely  exceeding  + .1.  The  quadrature  spectrum/co-spectrum 
ratio  is  correspondingly  small,  except  at  large  separations,  where  the 
co- spectrum  also  becomes  small.  Given  these  spectra,  the  assumption 
of  zero  net  vertical  energy  flux  appears  to  be  a justifiable  first 
approximation. 

The  Vertical  Squared  Coherence 

A normalization  of  the  vertical  cross  spectrum,  the  vertical 
squared  coherence 
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NORMALIZED  VERTICAL  CROSS  SPECTRUM 

HAWAII  NOV.  1973 
DEPTH  280  - 400  m 


5.63  11.25 


Frequency,  c.p.h. 

Figure  19.  Normalized  vertical  co-spectrum  (upper)  and  quadrature  spectrum  (lower) 
52  df.  Slashed  region  denotes  co-spectrum  between  . 7 and  . 8;  quadrature 
spectrum  between  .05  and  . 1;  dotted  region  denotes  co-spectrum  between 
.2  and  .3;  quadrature  spectrum  between  -.05  and  -.  1 
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VR2(az,Q)  = -v-^_ 

VC(o ,fi; 

is  most  frequently  used  to  gain  insight  into  the  distribution  of  en- 
ergy with  vertical  "mode".  Munk  and  Phillips  (1968)  using  sinusoids, 
have  identified  small  values  of  the  squared  coherence  in  a given 
frequency  band  with  the  relationship 

Akz  az  — 0(1) 

where  Akz  is  an  effective  bandwidth  of  energetic  wavenumbers.  Using 
Bessel  functions  as  representative  of  internal  vertical  wavefunctions , 
Garrett  and  Munk  (1972)  found  the  same  expression  generally  applicable. 

It  should  be  noted  that  the  assumption  of  random  phase  is 
necessary  in  the  development  of  the  above  relationship.  In  the  pre- 
sent case,  the  measurements  might  be  sufficiently  near  the  sea  surface 
that  the  phase  of  the  vertical  wavefunctions  might  not  yet  be  decor- 
related.  A more  precise  relationship  could  be  developed.  However, 
the  coherence  provides  information  primarily  about  the  major  spectral 
concentrations,  and  is  not  sensitive  to  secondary  concentrations  down 
many  db  from  the  main  peaks,  which  might  play  a di sproportional ly 
large  role  in  the  dynamics  of  the  wave  field.  As  a result,  major 

effort  will  be  expended  in  a more  direct  estimation  of  the  < depen- 
* 

dence  of  dE(r,0,ft). 

The  difficulty  with  the  coherence  technique  is  that  experi- 
mental noise  can  have  a strong  effect  on  the  measured  coherences, 


49 


Pinkel 


resulting  in  erroneous  bandwidth  interpretations.  Random  instrument 
noise  typically  acts  to  reduce  the  coherence  calculated  between  pairs 
of  probes.  In  the  present  case,  where  measurements  at  many  depths 
are  made  with  a single  instrument,  noise  tends  to  offset  the  isotherm 
displacement  estimates  uniformly,  increasing  the  measured  coherences. 
The  coherence  spectra  to  be  shown  display  significant  coherences  at 
frequencies  above  the  Vaisala  cutoff.  This  is  attributable  to  the 
dominance  of  thermometer  noise  in  the  data  in  this  frequency  range. 

At  frequencies  below  N(z),  where  there  is  a significant  internal  wave 
signal  (the  signal  to  noise  ratio  is  much  greater)  this  effect  is  not 
important.  The  error  effects  can  be  "reversed"  by  calculation  of  ver- 
tical coherences  between  isotherms  followed  by  separate  thermometers 
(Fig.  21).  Further,  the  "slant"  coherence  between  isotherms  followed 
from  separate  profilers  on  FLIP  was  calculated.  This  was  based  on 
information  taken  from  separate  vibratrons  and  therometers.  V ibratron 
noise  contributed  insignificantly  to  the  difference  between  the 
vertical  and  slant  coherence. 

Contour  plots  of  the  squared  coherence  estimate  versus  stret- 
ched depth  separation  and  frequency  are  given  in  Figure  20.  The 
actual  pattern  of  the  contours  is  quite  striking,  and  represents  a 
significant  departure  from  past  observations.  At  all  frequencies 
where  there  is  substantial  wave  energy,  the  coherence  is  above  .5 
(coherence  squared  above  .2b)  out  to  a separation  of  50  m.  The 
coherence  associated  with  a given  vertical  separation  tends  to  increase 
slowly  and  irregularly  with  increasing  frequency  to  2.5  cph.  Then  the 
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increase  becomes  extremely  rapid.  From  4 cph  to  the  Vaisala  cutoff 
the  vertical  motions  in  the  upper  thermocline  are  significantly 
coherent  for  depth  separations  in  excess  of  170  meters.  Above  N(z), 
the  coherence  drops  abruptly.  Its  behavior  is  determined  by  the  above 
mentioned  thermometer  noise  (Fig.  20)  or  the  statistical  instability 
of  the  estimate  (Fig.  21). 

The  high  values  of  the  squared  coherence  and  the  trend  toward 
increasing  coherence  with  increasing  frequency  (for  fixed  separation) 
sharply  contrast  with  much  of  the  vertical  coherence  information  ob- 
tained by  previous  workers  (Seidler,  White).  Particularly  disparate 
is  the  pioneering  study  of  Webster  (1972)  who  used  moored,  fixed  depth 
current  meters.  His  observations  indicated  that  the  depth  separation 
associated  with  the  .5  coherence  contour  should  decrease  with  increas- 
ing frequency,  reaching  a separation  of  13  meters  at  1 cph.  Garrett 
and  Munk  (1972)  have  attributed  this  discrepancy  to  the  microstructure 
contamination  of  Webster's  fixed  depth  instruments. 

Recent  results  of  Cairns,  measuring  temperature  from  a profil- 
ing midwater  float  (personal  communication)  and  Briscoe,  measuring 
temperature  from  a surface  spar  buoy  (personal  communication)  tend  to 
substantiate  the  general  pattern  seen  in  this  experiment. 

The  vertical  coherence  results  suggest  that  high  frequency 
vertical  motions  in  the  upper  thermocline  are  dominated  by  motion  of 
small  vertical  wavenumber.  Lower  frequency  motions  are  associated 
either  with  larger  vertical  wavenumbers , a wider  range  of  energetic 
wavenumbers,  or  some  combination  of  the  two  effects.  Next  a direct 
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resolution  of  the  vertical  spectrum  will  be  attempted  in  hopes  of 
further  examining  the  above  uncertainty. 


Jhe  Wavenumber  Frequency  Spectrum  from  Vertical  Separation 

The  objective  of  this  section  is  the  estimation  of  the  wave 
number  frequency  spectrum 

Mr 

dE(<»w)  = f dECce.J 

o=o 


from  the  depth/time  variability  of  the  wave  field.  The  motivation  for 
tins  effort  can  be  seen  in  Figure  22,  in  which  the  complex  isotherm 
displacement  Fourier  coefficients  n(*0,z,Q)  in  a given  frequency 
band  are  plotted  versus  depth.  The  smooth  undulation  of  the  Fourier 
coefficients  with  depth  is  the  inducement  to  attempt  a vertical  spec- 
tral analysis  of  the  data. 

The  procedure  to  be  used  is  considerably  more  complex  than 
standard  time  series  analysis  techniques  due  to  the  depth  variability 
of  the  Vaisala  frequency.  There  are  two  important  operational  con- 
sequences of  this  nonstationarity.  First,  unlike  standard  time  series 
analysis,  the  potential  energy  is  not  directly  given  by  the  total  in- 
tegrated variance  of  the  measured  signal.  Statistical  operations 
which  conserve  potential  energy,  not  variance,  must  be  employed. 

Second,  sinusoids  ate  not  the  natural  set  of  orthogonal  functions 
appropriate  tor  the  decomposition  of  the  signal.  Fourier  analysis 
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would  result  in  vertical  "wavenumber"  coefficients  which  are  not  sta- 
tistically independent.  Indeed,  the  concept  of  a vertical  wavenumber 
is  meaningful  only  in  the  W-K-B-sense.  In  the  analysis  outlined 
below,  the  necessity  for  dealing  with  vertical  wavenumber  will  be 
avoided;  equivalent  horizontal  wavenumber  magnitude 


k„ * kz<‘>  VH 

fi  - f 


2 -1/2 


will  be  utilized  instead. 

For  each  horizontal  position  x.,  there  are  averaged  values 

J 

of  nUj.z^fi)  in  forty  depth  intervals  zr  These  complex  non- 
stationary depth  series  are  to  be  expanded  as  a sum  of  yet  unspecified 
orthogonal  functions. 


40 


with 

an  = ^ Fn^  Y^z)  dz 


where  D is  the  depth  interval  spanned  by  the  array,  F is  the 

th  , n 

n function  in  the  orthogonal  set,  and  Y is  a weight  function  to 

be  determined. 

The  functions  {F^}  should  be  chosen  such  that  the  resulting 
separation  coefficients  fan)  are  statistically  independent,  i.e., 
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<Vj*>  = <lai^>  5 


1J 


Also,  in  order  that  the  average  potential  energy  be  given  by  the  sum 
of  the  squares  of  the  coefficients  {an}  it  is  necessary  that 


40 


<la„|Z>  - V2  pD"1  | n2(z)  <|n(z,Q)|2>  dz. 


It  can  be  demonstrated  (Appendix  II)  that  the  functions 

rn^zi^  " Gn^zi  )/fJ(zi ) ’ where  Gn  are  the  eigenfunctions  of  the  N 
weighted  vertical  co-spectral  matrix, 


[N(zi)  rp  VC  (zi  ,z.,a)  N(z  ■ )] 


are  the  desired  orthogonal  functions,  under  weighting  N2(zi ).  The 


resulting  separation  coefficients  {an)  are  independent.  The  variance 
<lail  > is  9iven  by  the  ith  eigenvalue  of  the  matrix  and  the  total 


variance 


40 


X 


i=l 


ai 


is  proportional  to  the  r.m.s.  potential  energy. 

The  eigenvalues  of  the  weighted  co-spcctral  matrix  are  thus  a 
complf  l.r  energy  conserving,  statistical  description  of  this  linear, 
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nonstationary  process.  However,  only  in  special  cases  will  the 
functions  (F^z.)}  resemble  the  dynamic  functions  (in  this  case 
t,le  W^,Q,zi))  v/llich  are  physically  meaningful  in  describing  the 
process.  It  is  necessary  to  relate  the  distribution  of  potential 

energy  among  the  functions  {Fn>  (as  given  b^  the  eigenvalu  of  the 
matrix)  to  dE(,.,J. 

The  relationship  between  the  eigenvalues  and  the  energy 
spectrum  is  given  by 

an  = j n Fn  II2  dz 

= j j X {|  W(K,«fz)  Fn(z)  N2(z )dz} 


= { dA  eitc  X Cn(K ,n) 

K 0 

where 
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2 

If  Cn  (k,q)  were  a delta-like  function,  strongly  peaked  at 

some  < = kn>  the  \ could  be  taken  as  a weighted  estimate  of  the 

spectrum  near  kn>  In  the  present  case  the  form  of  the  eigenfunctions 

is  strongly  influenced  by  the  vertical  sampling  interval  and  the 
2 

functions  {Cn  are  not  sharply  peaked. 

A 

It  is  reasonable  to  compose  an  estimate  E^(n)  of  dE(K,^) 
as  a linear  sum  of  the  eigenvalues: 

40 

Ek(a)  = “kn  xn 


M«)  = 


40 

L 

n=l 


‘kn  °n 


C/  U,n)dE (k,q) 


The  a^n  are  usually  chosen  to  make  the  spectral  window 
» -)  - £ a|<n  cn  (<  »fi) 

as  delta-like  as  possible.  This  is  done  by  minimizing  the  square  of 
the  window  outside  of  a desired  passband. 

A 

In  this  case,  where  it  was  assumed  that  the  spectrum  dF.^,^) 
would  tend  to  decrease  with  increasing  k,  the  square  of  the  weighted 
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window  K"nWk(K',fi)  was  minimized  outside  of  the  passband.  This 
had  the  effect  of  minimizing  the  error  in  the  spectral  estimate 


W 


k(K,n)  dE(<  ,n) 


W (<,Q)/Kn 


assuming  a K~n  form  for  dE. 

Trial  sections  of  data  were  repeatedly  analyzed  with  windows 
weighted  by  a factor  <"2,  k“1,5,  k"1  and  <°.  The  weighting  had 
only  a slight  effect  on  the  slopes  of  the  resulting  spectral  estimates 
Ek(w).  A k 1 weighting  was  finally  chosen  as  it  tended  to  produce 
a slightly  smoother  spectral  estimate. 

For  each  cruise,  sixteen  weighted  vertical  co-spectral  matrices 
were  calculated  from  the  sample  cross-spectrum  for  sixteen  frequency 
bands,  .23  cph  to  3.75  cph.  The  analysis  was  not  performed  at  higher 
frequencies  as  the  high  vertical  coherences  observed  above  4 cph  sug- 
gested that  the  first  mode  was  strongly  dominant  in  this  region.  The 
eigenfunctions  and  eigenvalues  for  these  40  x 40  matrices  were  then 
obtai ned . 

From  measured  representative  Vaisala  profiles  for  each  cruise, 
theoretical  vertical  wavcfunctions  W(>  ,n, z)  were  numerically  calcul- 
ated by  integrating  the  linear  internal  wave  equation  with  the  initial 
condition  U = 0 at  z = 0.  For  each  frequency,  the  equation  was 
integrated  for  120  evenly  spaced  values  of  k.  The  wavenumber  incre- 
ment A I was  chosen  such  that  the  maximum  wavenumber  r ~ 120  Ak 

max 
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corresponded  to  a wave  function  which  had  twenty  zero  crossings  from 
the  surface  down  to  440  m.  The  functions  were  normalized  such  that 


N2W2dz  = 1 . 


Using  the  previously  calculated  eigenfunctions,  a 40  by  120 
array  of  coefficients 

cn(*.n)  - CnK(n)  = | n2(z)  Fn(z)  WK(n,z)dz 


was  calculated  in  each  frequency  band. 

The  coefficients  a were  then  found  which  minimized 

wk2(K,fi)/K2  = wkK(n)/<2 

outside  of  the  desired  passband,  subject  to  the  constraint  W,  = 1 

k< 

in  the  passband. 

Forty  spectral  estimates  Ek(ft)  were  produced  in  each  fre- 
quency hand  for  values  of  k = 3 * l * Ak,  i = 1 , 40.  The  k corres- 
pond to  wavefunctions  making  from  1/2  through  20  oscillations  in  the 

top  *10  n.  With  window  resolution  set  at  6 Ak,  adjacent  estimates  were 
dependent . A contour  plot  of  a representative  window  function 
t!'kK ( ) is  presented  in  Figure  23. 
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Spectral  estimates  are  plotted  for  several  n in 

Figures  24-26.  At  all  frequencies  the  dominant  energy  is  associated 
with  vertical  scales  long  compared  tn  the  vertical  length  of  the  array. 
The  extent  to  which  this  long  wavelength  energy  contaminates  the  rest 
of  the  spectrum,  given  by  the  spectral  window,  is  plotted  as  a solid 
line. 

Three  spectral  regions  are  apparent.  The  spectral  maximum 
occurs  at  long  wavelength.  At  low  frequencies,  less  than  ~1.5  cph, 
the  array  resolution  is  sufficient  to  limit  this  peak  to  the  first 
twelve  modes.  Above  "3  cph  the  first  two  modes  are  clearly  dominant. 
This  pattern  is  compatible  with  the  vertical  coherence  information 
previously  presented. 

The  long  wavelength  maximum  drops  at  an  approximate 

•2 

to  k slope  down  to  a spectral  "shoulder".  This  shoulder  region 
has  a k®  to  slope.  Only  in  the  lower  frequency  bands  is  it 

clearly  above  the  spectral  noise  levels. 

At  higher  wavenumbers , the  spectral  intensity  in  each  fre- 
quency band  drops  precipitously  and  is  quickly  contaminated  by  the 
sidebands  of  the  long  wavelength  spectral  peak. 

This  information  can  be  alternately  portrayed  as  a contour  map 
of  potential  energy  versus  frequency  and  horizontal  wavenumber  (Figs. 
27-29).  For  the  November  1972  data,  CALCOFI  hydrographic  data  were 
available  from  top  Lo  bottom  of  the  ocean.  Coupled  with  FLIP-bourd 
density  estimates  in  Uie  top  400  in,  information  was  available  to 
generate  numerically  the  linear,  shear- free  internal  wave  dispersion 
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Figure  24. 


k H cycles/  meter 

f 2 2 

|l/2D  ^PN  W dz  • Ek  (S2)]  November  1972.  Frequency  bandwidth  0. 23  cph; 
ergs/cm3/cph  cpm. 
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kH  cycles/  meter 

Figure  26.  [1/2D  N2  W2dz  • E^fS)]  November  1973.  Frequency  bandwidth 
0. 23  eph;  wavenumber  bandwidth  0. 58x  lO^cpm;  ergs/cm3/cph  cpm 
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June  1973:  Log  potential  energy  versus  frequency  and  wavenumber.  Contour  interval  is  5 db  down 
from  10^  ergs/cm^/cph  cpm.  Dotted  region  denotes  first  spectral  band,  illustrating  variation  of 
wavenumber  bandwidth  with  frequency.  Slashed  areas  are  prohably  contaminated  by  leakage  of 
lower  wavenumber  energy. 
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relation,  which  is  also  plotted  on  Figure  27. 

It  is  apparent  from  these  spectra  that  the  number  of  excited 
internal  wave  modes  is  a strong  function  of  frequency.  Above  3 cph, 
the  first  mode  clearly  dominates  the  spectrum.  Spectral  levels  are 
undetectably  low  above  the  first  few  modes.  Although  the  greatest 
amount  of  energy  is  in  the  lower  modes  (longer  horizontal  wave- 
lengths) at  lower  frequencies  also,  there  is  significant  short  wave- 
length energy,  associated  with  the  spectral  shoulder.  It  would  be 
difficult  to  infer  the  presence  of  this  short  wavelength  energy  from 
the  vertical  coherence  information  alone,  as  it  is  down  10-15  db  from 
the  long  wavelength  spectral  maximum. 

This  information  is  al ternatively  contoured  versus  frequency 
and  W.K.B.  equivalent  vertical  wavenumber 


where 


, 2 

k = k (—  - ~ J \ 
z KH  1 2 2' 

Si  - p 


V2 


| N2(z)dz 
0 


in  Figures  30-32.  The  frequency  dependence  of  F,  (o)  can  be  removed 
by  multiplying  by  u . When  this  is  done,  the  contour  lines  run 
nearly  parallel  to  the  o axis.  This  suggests  a vertical  wavenumber 
bandwidth  which  is  independent  of  frequency.  This  simple  pattern  is 
deformed  at  the  lnwc*rl  frequency,  corresponding  *r>  the  increased 
energy  in  the  spectral  shoulder.  The  agreement  with  the  Garrett -Munk 
1974  spectrum  is  good,  except  in  the  renion  of  the  spectral  shoulder. 
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Vertical  Wavenumber  cpm 


Figure  30.  November  1972:  Log  potential  energy  versus  frequency  and  W.  K.  B.  vertical 
wavenumber.  Contour  interval  is  5db  down  from  .32x  104  ergs/cm3/cph  cpn- 
Dotted  line  represents  10  mode  dispersion  line.  Slashed  region  is  probably 
contaminated  by  leakage  of  low  wavenumber  energy. 
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Vertical  Wavenumber  cpm 

Figure  32.  November  1973:  Log  potential  energy  versus  frequency  and  W.K.B.  vertical 
wavenumber.  Contour  interval  is  5 db  down  from  If)4  ergs/cm^/cph  cpm. 
Slashed  region  is  probably  contaminated  by  leakage  of  low  wavenumber  energy 
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Above  2 cph,  a fixed  vertical  wavenumber  bandwidth  corresponds  to  a 
d^cjgdsjng  number  of  inodes  with  increasing  frequency.  Th^s  is  re- 
flected in  the  increase  in  vertical  coherence  mentioned  in  the 
proceeding  section. 

It  should  be  pointed  out  that  energy  in  spectral  regions  of 
low  phase  velocity 

Cph  10  cm/sec 

KH 

might  be  significantly  Doppler  shifted.  The  energy  must  be  severely 
shifted  before  any  major  change  in  the  general  spectral  pattern 
occurs . 

It  is  important  to  check  these  spectral  estimates  for  consist- 
ency with  other,  independent  experimental  evidence.  Simulated  "towed 
spectral  estimates"  can  be  generated  from  Figures  27  and  29  by  inte- 
grating the  energy  in  each  band  and  dividing  by  kR  to  convert 
from  radial  to  linear  co-ordinates.  The  assumption  of  horizontal 
directional  isotropy  is  necessary  here.  The  simulated  towed  spectra 
arc  presented  in  figure  33.  They  have  a JkJ"2  slope  at  low  wave- 
number, which  steepens  slightly  for  wavenumbers  larger  than  .5  cycles 
per  kilometer.  These  simulated  results  compare  remarkably  well  with 
1 lu  towed  spectra  of  El i Katz  (personal  communication). 

A second  check  can  be  made  using  the  isotherm  slope  spectrum 
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z j j (k  cosq )2dE(K  ,q ,q)  W2(k,q,z). 

K 0 

Such  a spectrum  is  obtainable  from  measurement  of  the  tilt  of  isother- 
mal surfaces,  using  the  three  profiling  units  on  FLIP.  They  will  be 

p 

presented  in  the  next  section.  Due  to  the  K weighting,  this 
spectrum  is  very  sensitive  to  the  short  wavelength  energy  levels  in  a 
given  frequency  band. 

Figure  34  presents  the  simulated  slope  spectra  generated  by 
performing  the  weighted  integration  of  E(. (n).  The  spectral  levels 
are  underestimates,  as  the  k integration  was  terminated  at  large  k 
when  the  vertical  spectral  window  indicated  that  the  quality 

of  E^.(j)  was  suspect.  This  potential  underestimation  is  most 
significant  in  the  lowest  and  highest  frequency  bands.  Horizontal 
isotropy  was  again  assumed,  to  perform  this  conversion. 

An  approximate  f>  slope  is  given  by  the  simulation.  This 
contrasts  with  the  Garrett  and  flunk  (1972)  spectrum  which  predicts  a 
white  slope  spectrum,  independent  of  the  number  of  modes  chosen. 
Similarity  between  the  slope  spectral  estimates  to  be  presented  and 
the  simulated  spectra  will  give  confidence  in  the  accuracy  of  the 
vertical  separation. 

The  I s o the  nil  Slo  p o Spec  t rum 

The  isotherm  slope  spectrum 


w 
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SLx(z,ft)  = <||^  (x,z,n)  | 


i s 


given  in  terms  of  the  equivalent  continuum  representation  by 


SLx{z,n)  j j ( k cose)  <dAdA*  (K,o.n)>  Vt  {Kyntz) 
K o 


An  approximation  to  the  slope  spectrum  can  be  obtained  from  the 
present  data  by  assuming 


gG'Z.a)  ~ 1 n(X2’2,n) 

IX1  ' X2  I 


ax 


The  approximate  slope  spectrum  has  the  form 


SLx(z,o)  = T-T  = 1 


f 2,  ,1*1  ‘ X2 


(1  - *2r  i > 


sin  (<(——= cose)) 


<dAdA*(k,o,fi)>  W%:,n,z) 


For  wavelengths  long  compared  to  the  probe  separation 


. 2/  IX1  " X?  ? X1  X9  ^ 9 

sin  (k  2 cosq)~  k 2J cos  o 


and  the  approximation  is  good.  Were  appreciable  energy  concentrated 
in  wavelengths  of  the  same  order  as  the  array  separation  lennth  or 
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smaller,  the  approximate  slope  spectrum  would  tend  to  underestimate 
the  true  spectrum. 

Typical  approximate  slope  spectra  are  presented  from  the 

November  1972  cruise  in  Figure  35.  The  slopes  are  taken  with  respect 

to  the  mean  downwind  and  mean  crosswind  directions.  Below  1 cph  the 
-1/2 

spectrum  has  an  q frequency  dependence.  Above  this  frequency 
there  is  a spectral  shoulder  of  lesser  slope.  At  frequencies  far 
above  the  Vaisala  frequency,  the  spectrum  is  white.  A pronounced 
cutoff  is  not  apparent.  There  is  no  apparent  systematic  variation  of 
the  slope  spectrum  either  with  direction  or  depth  which  is  present 
in  all  three  cruises. 

The  general  form  of  the  slope  spectrum  is  quite  similar  to 
the  simulated  slope  spectrum  inferred  in  the  previous  section  from 
the  vertical  variabilly  of  the  data.  The  discrepancy  between  the 
spectrum  and  the  Garrett-Munk  1972-74  models,  which  predict  an  n° 
spectral  form  might  reflect  the  slight  increase  in  the  low  frequency, 
high  mode  energy  discussed  in  the  previous  section,  and  not  included 
in  the  models.  Garrett  and  Munk  emphasize  that  this  is  an  "isotherm 
depth  difference"  spectrum,  not  a true  slope  spectrum.  The  Garrett- 
Munk  1974  isotherm  difference  spectrum  has  an  approximate  slope 

in  the  intermediate  frequency  range  and  agrees  reasonably  well  with 
this  spectrum  (Garrett  and  Munk,  personal  communication). 

Aj i r.stimate  j >c  Hnr iron t ,:1  l.'avcnutnhcr 

fi'iure  16  presents  an  isotherm  displacement  spectrum  and  an  isotherm 
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difference  spectrum  from  the  November  1972  cruise.  A comparison 
of  the  spectral  forms  provides  some  information  on  the  wave- 
number  dependence  of  the  spectrum. 

A horizontal  wavenumber  estimate  can  be  defined 


<x(z,fi)  = [SLx(z,n)/S(z,Q)]1/2 


(kcos^)2  <dAdA*>  W2 


1/2 


K 0 


f f 2 

1 I <dAdA*>  VT 


K 0 


for  each  frequency  and  depth.  The  relationship  between  this  horizon- 
tal wavenumber  estimate  and  the  wavenumber  dependence  of  the  spectrum 
is  precise  only  in  the  case  of  uni-modal  propagation 


dE(r,0,(,i)  = d E ( 6 , (o ) 6 {y  - <(co)). 


Here  r = k(w). 

It  is  useful  to  consider  two  other  cases.  Were  the  spectrum  in 
a given  frequency  band  directionally  isotropic,  constant  between 
wavenumbers  <(|^n  and  <max>  and  zero  elsewhere,  the  resulting  scale 
wavenumber  magnitude  would  be. 


- (r  2 + r 2)1/2  = (1  (k  2 4 ■ K . y + . 2 ) ) 1/2 

> y 3 max  nnn  max  ‘nnn  '' 
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A second  isotropic  spectral  model,  with  a k dependence  between 
wavenumbers  K^n  and  <r|ax  and  zero  elsewhere,  is  associated  with 
scale  horizontal  wavenumber 

< = (k  r 11//2 

kH  umin  "max' 

Figures  37-39  portray  versus  frequency  for  the  November 
1972  cruise.  Isotherm  information  from  mean  depth  220  m was  used. 

The  first  ten  dispersion  lines  are  also  drawn  on  the  plots.  Below 
3 cph  the  wavenumber  increases  slowly  with  frequency,  then  more 
rapidly,  following  the  first  mode  dispersion  line  When  n exceeds 
N,  the  estimate  becomes  erratic. 

Using  the  models  introduced  above,  an  estimate  can  be  made 
of  the  bandwidth  of  energetic  wavenumbers.  If  it  is  assumed  that 
Kmin  is  given  by  the  first  mode  dispersion  line,  then  the  variation 
°f  Kmax  1S  fetched  in  Figure  40.  The  dominance  of  the  first  mode 
above  3 cph  and  the  general  widening  of  the  wavenumber  bandwidth  be- 
low are  apparent  in  both  models.  However,  due  to  the  oversimplified 
form  of  the  assumed  wavenumber  dependencies,  neither  model  depicts 
the  wavenumber  bandwidth  variation  accurately  in  the  lower  frequency 
region. 

Ho rizontal  Coherence 

Horizontal  coherence  measurements  have  been  utilized  by 
other  experimenters  to  provide  an  indirect  estimate  of  the 
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f 


Figure  38.  November  1972:  Crosswind  component  of  horizontal  wavenumber 
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horizontal  wavenumber  dependence  of  the  spectrum.  In  the 
present  case,  these  measurements  are  of  somewhat  limited  value.  The 
horizontal  dimensions  of  the  array  are  sufficiently  small  compared 
to  typical  internal  wave  scales  that  high  coherences  are  measured  at 
all  frequencies  where  waves  are  present.  The  significant  information 
on  the  wavefield  is  given  by  small  deviations  of  the  coherence  from 
1.  A noise/signal  ratio  which  is  acceptable  for  measurements  where 
lower  coherences  are  expected  might  strongly  affect  the  interpreta- 
tion of  high  coherence  measurements. 

Figure  41  presents  an  estimate  of  the  horizontal  coherence 
squared  between  5.6  hour  records  of  the  14°  isotherm  depth  fluctua- 
tion as  followed  by  the  port  and  starboard  probes  (separation  44  m). 
This  particular  data  run  was  rapidly  sampled,  the  Nyquist  frequency 
being  45  cycles  per  hour.  The  associated  displacement  spectrum  is 
given  in  Figure  8.  The  squared  coherence  is  statistically  insignifi- 
cant above  the  Vaisala  frequency  and  nearly  unity  from  the  cutoff 
down  to  approximately  three  cycles  per  hour.  Below  this  frequency 
the  value  drops  corresponding  to  the  increased  horizontal  wavenumber 
bandwidth  discussed  in  the  preceeding  section.  Presumably  the 
squared  coherence  should  then  rise  toward  lower  frequencies.  The 
poor  resolution  and  stability  of  this  record  obscures  any  signifi- 
cant low  frequency  rise. 

A contour  plot,  of  horizontal  coherence  squared  versus  depth 
and  frequency  is  presented  in  Figure  42.  Three  days  of  data  from  the 
Hawaii  cruise  were  utilized  to  produce  this  ?P,  degree  of  freedom 
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Frequency,  c.p.h. 


STAR  BOA  11  D & AFT  PROBES 


Frequency,  c p.h 

Figure  42.  Horizontal  coherence  squared  versus  depth  and  frequency,  28  df. 

Shaded  region  represents  coherence  squared  between  . 7 and  . 8. 
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estimate.  The  depth  variability  of  the  high  frequency  Vaisala  cut- 
off is  quite  apparent.  The  squared  coherence  between  the  starboard 
and  aft  probes,  which  lie  along  FLIP'S  drift  path,  tend  to  cut  off  at 
a slightly  higher  frequency  than  between  the  starboard  and  port 
probes,  which  are  somewhat  more  perpendicular  to  the  drift  direction. 
This  effect  could  be  a manifestation  of  the  Doppler  shifting  of  the 
measurements,  or  a result  of  short  term  (3  day)  anisotropy  of  the 
wave  field. 

Given  the  marginal  statistical  stability  of  these  estimates, 
only  the  most  general  conclusions  .an  be  drawn.  Below  3 eph,  the 
squared  coherence  across  the  array  is  approximately  .7-. 8,  indicating 
a relatively  wide  range  of  horizontal  wavenumbers.  Above  3 eph  the 
squared  cohererce  tends  to  be  greater  than  .9,  suggesting  a very 
narrow  range  of  energetic  wavenumbers.  This  supports  the  vertica1 
coherence  measurements,  which  strongly  suggest  that  many  modes  are 
important  at  the  low  frequencies,  while  the  first  few  modes  dominate 
the  spectrum  in  the  higher  frequency  region. 

An.  Invest i_g a_ti on  o f_Jn ternal  Wave  D i spe rs ion 

In  the  discussion  of  the  proceeding  sections,  it  has  been 
assumed  that  the  observed  motions  were  linear  internal  waves.  An 
attempt  will  now  be  made  to  establish  the  validity  of  this  assumption. 

It  should  be  noted  that  the  general  analysis  approach  has  not 
been  dependent  on  the  fact  that  the  motions  are  internal  waves.  The 
vertical  separation  parameter  k^  was  introduced  as  an  independent 
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variable  characterizing  the  oscillations  of  a convenient  set  of 
orthogonal  functions,  Wk(u,z).  Only  if  linear  internal  wave  theory 
well  relates  the  horizontal,  vertical  and  temporal  scales  of  the 
motions,  should  the  parameter  kH  physically  correspond  to  a hori- 
zontal wavenumber.  Using  information  from  the  three  profilers  on 
FLIP,  direct  estimates  of  horizontal  wavenumber,  r,  will  now  be 
calculated,  for  each  fl.k  band.  In  those  regions  of  the  wavenumber 
frequency  domain  in  which  i - kH.  linear  internal  wave  theory  pre- 
sumably provides  a meaningful  approximate  characterization  of  the 
motions. 

Central  to  this  discussion  is  the  vertically  separated  hori- 
zontal cross  spectrum 


which  is  obtained  in  a manner  analogous  to  the  vertically  separated 
wavenumber  spectrum  Ek(n)  (Appendix  III).  Combining  these  two 
spectra,  a vertically  separated  slope  spectrum 

SLk(“’Axij)  = <2  Ek(n)  - 2 Real  (VCRjfi.Ax. .)))  ax..  "2 

^ J 1 J 

can  be  obtained.  The  ratio  of  the  slope  spectrum  to  the  energy 
spectrum  will  be  used  to  estimate  a horizontal  wavenumber,  r,  in 
the  manner  discussed  previously. 


92 


SIO  Reference  75-10 


Horizontal  wavenumber  magnitude  estimates  k( k ,n)  are 
presented  for  several  frequency  bands  in  Figures  43-45 . The  vertical 
scale  in  the  figures  is  not  kH  but  the  W.K.B.  mean  vertical  wave- 
number 


N2  - v2?'2 
kz  = kH  ^2  _ ^.2 


Reference  lines  K = kH  are  also  plotted  for  each  frequency  band. 
Agreement  between  measured  wavenumbers  and  linear  theory  is  quite 
good.  The  wavenumber  ratio  kz/<  is  approximately  constant  at  a 
given  frequency,  as  emphasized  by  Cckart  (I960).  The  frequency 
variation  of  the  wavenumber  ratio 


r,2  2 

kz'  * ' 2 

0)  - f 


1/2 


is  also  clearly  apparent. 

The  horizontal  beamwidth  of  the  wavefield  can  be  related  to 
the  vertically  separated  squared  horizontal  coherence 


Rk  <"•%> 


|VCRk  (il.ix.j)!2 

Tk  («.*,)  tknrx-) 


as  discussed  by  Munk  and  Phillips  (1968).  For  small  k,  the  co- 
herence across  the  array  was  found  to  be  quite  high  at  all 
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Figure  44.  June  1973:  Horizontal  wavenumber  magnitude  versus  W.K.  B.  vertical 
wavenumber.  Solid  line  gives  theoretical  relation  ku=x. 
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frequencies.  With  decreasing  vertical  wavelength  it  dropped  rapidly. 
The  coherence  was  frequently  below  that  expected  from  an  isotropic 
spectrum  with  energy  concentrated  at  wavenumber  k^.  The  values  were 
consistent  with  a modeler  isotropic  spectrum  with  energy  restricted 
to  some  thin  annular  nng  about  k^.  This  is  reasonable,  considering 
the  finite  resolution  of  the  vertical  array.  Instrument  noise,  im- 
precision in  the  analysis,  and  other  experimental  difficulties  must 
also  play  a role  in  reducing  the  coherences. 

Given  the  low  coherences,  there  is  no  way  to  determine  any 
departure  from  isotropy  in  the  signal.  Therefore,  no  extensive  effort 
was  made  to  determine  the  0 dependence  of  the  spectrum. 

A previous  study  of  the  high  frequency  waves  (ft  > 3 eph), 
not  presented  here,  demonstrated  significant  anisotropy  over  1 day 
time  periods.  Directional  beamwidths  of  90°  to  180°  were  found  to 
cnaracterize  the  predominantly  first  mode  motions. 
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CHAPTER  IV 
CONCLUSIONS 

Above  a minimum  horizontal  phase  velocity,  5-10  cm/sec,  the 
observed  isothermal  motions  could  be  well  represented  by  a continuous 
spectrum  of  internal  waves.  The  dispersion  properties  of  these 
motions  agreed  with  linear  theory  to  within  experimental  accuracy. 
Below  this  minimum  phase  velocity,  which  corresponds  to  a minimum 
vertical  wavelength  *z  40  in  at  low  frequency,  Az  ^100  m at  high 
frequency,  the  motions  could  not  generally  be  successfully  decomposed 
in  terms  of  internal  wave  vertical  wave  functions  exhibiting  the 
dispersion  properties  predicted  by  linear  theory. 

At  frequencies  far  above  inertial,  the  internal  wave  spectrum 
can  be  divided  into  two  regions,  according  to  frequency.  Horizontal 
coherence  measurements  suggest  that  the  horizontal  wavenumber  band- 
width is  narrow  above  *'2  eph,  somewhat  wider  below.  In  the  lower 
frequency  region,  the  frequency  dependence  of  the  spectrum  is  q"^. 
Both  the  vertical  coherence  and  horizontal  slope  spectrum  suggest 
that  the  wavenumber  bandwidth  decreases  with  decreasing  frequency. 

In  this  frequency  range,  the  horizontal  wavenumber  dependence 
oi  fhe  wavenumber  f requeney  spectrum  c«ni  be  discussed  in  terms  of 
three  distinct  regions.  At  the  lowest  wavenumbers,  there  is  a 
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slope.  This  levels  into  an  approximate  - k-1  spectral  shoulder. 
The  spectrum  then  decreases  rapidly  again  at  higher  wavenumber  into  a 
region  that  is  dominated  by  analysis  noise. 

These  measurements  do  not  preclude  the  Garrett  and  Munk  (1972) 
suggestion  that  there  is  a fourth  region  of  slope,  which 

embraces  the  first  several  modes  (which  cannot  be  resolved  by  this 
array). 

The  discovery  of  the  region  of  short  wavelength, 

low  frequency  energy  is  one  of  the  more  important  results  of  this 
experiment.  The  prominence  of  the  "shoulder"  appears  to  decrease 
with  increasing  energy  in  the  wavefield.  This  spectral  adjustment 
possibly  acts  in  such  a manner  that  equilibrium  shear  levels  in  each' 
frequency  band  are  maintained.  Low  frequency,  high  mode  energy  has 
been  suggested  by  Kenyon  (1968)  to  result  from  internal  wave  nonlinear 
interaction.  Cox  and  Sandstrom  (1962)  have  predicted  the  existence  of 
very  high  mode  energy  as  a result  of  the  sea  floor  scattering  of 
barotropic  motions.  On  the  other  hand,  this  shoulder  might  be  the 
spectral  signature  of  horizontally  advecting  features  in  the  temper- 
ature field  which  do  not  correspond  to  density  perturbations.  How- 
ever, as  the  motions  marginally  satisfy  the  internal  wave  dispersion 
relation  in  this  spectral  region,  it  is  doubtful  that  advective 
effects  have  completely  dominated  the  spectrum.  Further  work  will 
have  to  be  done  to  understand  the  role  of  this  low  frequency,  short 
wavelength  energy  in  the  dynamics  of  the  sea. 

Above  ~2  cph,  the  waveguide  thickness  decreases  with 
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increasing  frequency.  In  this  region,  the  first  mode  becomes 
increasingly  prominent  with  greater  frequency,  as  reflected  in  the 
vertical  coherence  calculations  as  well  as  the  vertical  spectral 
separation.  The  displacement  spectrum  and  slope  spectrum  have 
irregular  slopes,  due  to  the  intermode  interference  discussed  in 
Chapter  III. 

Given  the  low  vaiues  of  horizontal  oherence  measured  across 
the  array  and  the  long  horizontal  wavelengtis  of  the  motions  studied, 
it  is  not  possible  to  detect  any  significant  degree  of  directionality 
in  the  waveficld  over  10  day  periods.  However , this  conclusion  of 
isotropy  might  result  from  experimental  noise  which  also  acts  to 
reduce  the  coherence. 

Although  not  discussed  in  this  work,  significant  anisotropy 
was  detected  in  the  high  frequency  (u  > 3 cph)  motions  over  1 day 
records. 

A further  experiment  is  suggested,  utilizing  measurements  of 
density,  not  temperature,  and  sufficient  flowmeter  lowerings  such 
that  a clear  picture  of  the  mean  current  shear  can  be  formed.  The 
study  of  relationships  between  slope  spectral  levels,  mean  shear,  and 
regions  in  which  the  density  profile  is  rapidly  evolving  would 
greatly  increase  understanding  of  upper  ocean  dynamics. 
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x 

2 

t 

u(x,z,t) 

w(x,z,t) 

n(x,z,t) 

D 

H 

T(x,z,t) 

p(z) 

N(z) 

c 

CO 

n 

f 

K 

CP 

W(k,u,z) 

k 


kH  = a(u)k 


LIST  OF  SYMBOLS 

horizontal  coordinate 
vertical  coordinate 
time 

horizontal  velocity 
vertical  velocity 
vertical  displacement 

depth  interval  spanned  by  array  - typically 
60-440  m. 

depth  of  ocean 

temperature 

density 

Vaisala  frequency 

velocity  of  sound 

intrinsic  frequency 

encounter  frequency 

inertial  frequency 

horizontal  wavenumber 

horizontal  phase  velocity 

linear,  shear  free  internal  wave  vertical 
wavefunction 

dimensionless  independent  parameter  characteriz- 
ing the  oscillatory  behavior  of  a set  of  vertical 
internal  wavefunctions  Wk(u>,z)  used  to  decompose 
the  vertical  variability  of  the  motions 

dimensional  equivalent  horizontal  wavenumber 
magnitude:  Wk(w,z)  = W(kH,w,z) 
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SL-)  '/2 

9/ 


N2 

Kz  72 

0 - f 

ep(z) 

Ek!z) 

C (x-j  ~ X£»  z-j  ,Z2»  t-j 
CR(Xi  ~ X£»  z1  ^2*^) 

CR(x  - %2*  .z^,  n) 

dE(*;,G,w) 

2tt 

dE(K,u)  = ( dE(K,o,u) 


Ek(fi)  cph"1  cpm"1 

wk(K,fi)  = wkK(o) 


W.K.B.  equivalent  vertical  wavenumber 

r.m.s.  potential  energy 

r.m.s.  kinetic  energy 

t2)  process  covariance  function 

process  cross  spectrum 

sample  cross  spectrum 

three  dimensional  internal  wave  spectral 
function 


two  dimensional  internal 
wave  spectral  function 

estimate  of  dE(*,u)  in  frequency  band 
a,  vertical  separation  band  k 


vertical  spectral  window 
S(z,ft)  = CR(0,z,z,fi)  m^/cph  vertical  displacement  frequency  spectrum 
V(z,ft)  = O2  CR(0,z,z,o)  (cm/sec)2/cph 

vertical  velocity  frequency  spectrum 
vertical  cross  spectrum 


VC(z1,z2,fi)  = CR(0,z1 ,z2,n) 
VC(az ,n) 


vertical  cross  spectrum  vs  frequency  and 
W.K.B.  stretched  separation 
z< 


AZ  = - 
N 


1 z,2 


N(z)dz 


VR2(az,o)  = 

lycfeox 

VC(o,o)2 

vertical  squared  coherence 

SL-  (z.fi) 

cph"1 

Isotherm  slope  spectrum 

SH-  (z.n) 

-2  -1 
sec  cph 

horizontal  shear  spectrum 
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N 


VCRk  (n.AXjj) 


vertically  separated  horizontal  cross 
spectrum 


SLk(sj,Ax.j) 

Rk(s!-4xij) 


vertically  separated  slope  spectrum 

vertically  separated  horizontal  squared 
coherence 
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APPENDIX  I 

THF  STATISTICAL  DESCRIPTION  OF  THE  MOTION 


The  data  set  can  be  uniquely  described  by  the  Fourier- 
Stieltjes  representation 

to  2tT  » 

n(3T,7Jt)  = / / { dA(K,o,ui,z)e1'  K'x  " 

K=0  0 =0  U)=to 

where  n is  the  vertical  displacement  of  an  isothermal  surface  about 
its  mean  position  and 

dA(»c  ,o  ,z)  = dA(<  ,e  + v , -w,z) 


in  order  that  n be  real  valued. 

The  isotherm  displacement  variance  is  given  by 


<n2(z)' 


► 0 w 


<dA(<  ,f)  ,z)dA(£  ,tt  ,z)> 


ei  (<  + t)-x  - (u  + n)t. 


If 


,2(z)> 


is  to  be  stationary  in  x and  t. 


it  is  necessary  that 


1-1 


V « ^ 


^ I in  I I .I  f Mil.  W 1.  I IN.  «■  0« ’ ■"***'»  Ni>,l  Mr 


Pinkel 


^ ~ /W 

<dA(K,o,w,z)dA(«,,<{),n,z)>  = <dA(K,0,q,z)dA*(K,e,q,z)> 

6(k  + £)d( 0)  + ft)  . 

The  motions  of  primary  interest  to  this  study  are  those  with 
temporal  frequency  u less  than  the  local  Vaisala  frequency  N(z). 

In  this  case,  it  is  reasonable  to  separate  the  vertical  dependence  of 

— # 

dA( k, 0 , u), z ) = dA(<, g,oj)  W(<,(j,z) 

where  W is  a solution  of  the  linear,  shear  free  internal  wave 
equation 


i*  + k2  uIm  - -2)  w 

3?  (>  - f2) 


= 0 


where  W = 0 at  z = 0 and  is  normalized  such  that 


J 

2lJ 


p(z)  N2{z)  W2(K,u,z)dz  = 1 


Now  the  variance  is  given  by 


<n  (z)>  = 


<|dA|?,  U2(,,,u,z) 


J i ) 

k 0 o; 
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With  the  additional  r striction  that  W = 0 at  a flat  bottom, 
z = - II,  the  energy  is  quantized  into  discrete  modes 


<n2(z)>  = l 
i 


<|dA . (o,to)  |2>  W2(k.  ,to,z) 

0 co  1 


where,  for  each  frequency,  the  {W(«. ,w,z)}  are  an  orthogonal  set  over 

. N2  - co2 

the  interval  (0,  - II)  under  the  weight  function  (-^ . 

to  - f 

Garrett  and  Hunk  (1972)  have  suggested  that  the  effects  of  the 
particular  boundary  conditions  at  the  sea  surface  and  bottom  are  gen- 
ially unimportant.  They  then  introduced  an  "equivalent  continuum" 
representation 


n(x,z  ,t) 


j J dA(k-,0,<o)  W(<,to,z)e^  K x 

r 0 to 


tot 


such  that 

j „ 

< | dA j 2>  W2(k,,j,z)  = J.  < | dAi  1 2>  W2(<i,  to, 7.)  . 

\r  . 

i i 
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APPENDIX  II 

CHOOSING  THE  ORTHOGONAL  FUNCTION  SET 
For  each  horizontal  position  x.  there  are  measured  values 

J 

of  n(x.,z.,ft)  at  forty  depths  z..  This  complex,  ncnstationary 
J ^ 

depth  series  is  to  be  expanded  as  a sum  of  orthogonal  functions 

40 


wi  th 


a 


J n(z)  Fn(z)  Y(z)dz 
D 


where  F is  the  nth  function  in  the  desired  orthogonal  set  and 
n 

Y is  a weight  function  to  be  determined. 

In  order  that  the  potential  energy  be  given  by  the  sum  of  the 
squares  of  the  coefficients  {an),  ^ necessary  that 


40 

j,  <Un|2>  - 1/2  p D-1 


N2  <1 


>dz 


- 1/2  p D"1 


a * F -dz 

m m 


n-i 
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= 1/2  p D 


m 


dz  . 


2 

This  requirement  can  be  satisfied  if  <a„  a *>  = <|aj  > <$  m 
M £ m 1 m1  jtm 


or  N2  F F dz  = 6 where  the  weight  function  Y is  seen  to 
1 i m i,m  3 

D 

be  N2. 

It  will  now  be  demonstrated  that,  if  {Gn(z^)}  is  the  set  of 
eigenfunctions  of  the  N weighted  co-spectral  matrix 


<N(zi)  n (xQ,z .,n)  n*(x0,zj.,n)  N(z^)> 
then  the  functions 


<Fn<2i»  ' <Gn(z1)/N(zi)} 

are  the  desired  complete  orthogonal  set. 
Wri  te 


N(z,)  n(x0,Zj.Q)  = l an  Gn(Zj) 


a 


n 


j N(z.)  n(zi)  Gn(z.)dz  • 
D 
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The  N weighted  co-spectral  matrix  can  be  written 


<N(z.) 


W Gn<zj> 


where  from  linear  algebra  we  have 


l Gn(zi!  G*  (zi>  = 5 


r\z 


and 


<an  " "Mnl  ” 6 


nj i 


An5n* 


where  > is  the  nth  eigenvalue  of  the  matrix.  By  substitution 
F(z)  - G(z)/N(z),  it  is  seen  that 


^(xj*2i*“)  = l an  Fn(zi) 


an  = 


N2  n Fn  dz 


and 


D 


F F trdz  = fi 
n i R£ 


If  the  weighted  co-spectral  matrix  has  no  zero  eigenvalues, 
the  { G (z  )}  are  a complete  set.  It  can  easily  be  verified  that 
the  Fn(z.)  are  also  complete,  provided  that  N(z)  > 0 for  all  ; 
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APPENDIX  III 

A PRIORI  OPTIMAL  CROSS  SPECTRAL  SEPARATION 


Model  the  signal  as  a sum  of  "modes" 


120 


and  represent  it  as  a slum  of  eigenfunctions 

40 

ntZj.Xj)  ■ an(xj)  fn(2j) 


where 

an^xj ) = { n(zi«xj)  fn(zi)  Weight  (x. )dz 
D 


( 1?0 

= J JL  Am  (xj}  Fm(zi}  Vzi>  Weight(zi  )sz 


120 


= 1 
m=l 


A (x,)C 
mv  J mn- 
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and 

f,nn  = 1 Fm  fn  Wel'3ht  dz  • 
D 


Now 


120  120 


C f 
mn  Tin 


The  quantities 


can  be  calculated  from  which  we  wish  to  estimate 


< A. 


m 


(X,) 


A * 

m 


The  estimate  is  formed  as  a linear  sum 
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40 

= l 
£=1 

V 

40 

120 

= l 
1=1 

°W 

J,  <V*i>  V (xj)> 

120 

= l 

n=l 

W 

mn 

<Anfxi)  <*j)> 

The  cross  spectral  window  is 
40 


Wnn-  l 


1=1  amt  nt 


which  is  the  same  window  that  is  associated  with  the  vertical  spectral 
separation  in  Chapter  III.  Thus  the  same  coefficients  v which 

were  used  for  the  optimal  vertical  separation  can  be  used  to  form  the 
optimal  cross  spectral  estimate. 
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